Fluxon-semifluxon interaction in an annular long Josephson 0-7r-junction 



E. Goldobin, 1 ^ N. Stefanakis, 2 D. Koelle, 1 and R. Kleiner 1 

1 Physikalisches Institut - Experimentalphysik II, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 
2 Institut fur Theoretische Physik, Universitat Tubingen, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 
(Dated: February 2, 2008) 

We investigate theoretically the interaction between integer and half-integer Josephson vortices 
(fluxons and semifluxons) in an annular Josephson junction. Semifluxons usually appear at the 
0-7r-boundary where there is a 7r-discontinuity of the Josephson phase. We study the simplest, but 
the most interesting case of one 7r-discontinuity in a loop, which can be created only artificially. We 
show that measuring the current-voltage characteristic after injection of an integer fluxon, one can 
determine the polarity of a semifluxon. Depending on the relative polarity of fluxon and semifluxon 
the static configuration may be stable or unstable, but in the dynamic state both configurations 
are stable. We also calculate the depinning current of N fluxons pinned by an arbitrary fractional 
vortex. 

PACS numbers: 74.50.+r, 85.25.Cp 74.20.Rp 

Keywords: Long annular Josephson junction, sine-Gordon, half-integer flux quantum, semifluxon, 0-pi- 
j unction 



I. INTRODUCTION 

For conventional Josephson junctions the first Joseph- 
son relation reads I s = 7 c sin(0), where / s is the su- 
percurrent through the junction, I c is the critical cur- 
rent and cf> is the so-called Josephson phase which is 
equal to the difference of the phases of the quantum- 
mechanical macroscopic wave functions in the electrodes. 
The Josephson relation for a Josephson 7r-junction is 
I s = —I c sin((/>) = I c sin(0 + 7r), i.e., a 7r-junction can be 
considered as a junction with negative critical current or 
having an additional phase shift of 7r between the phases 
of the wave functions (therefore the name). Accordingly, 
conventional Josephson junctions are sometimes called 
0-j unctions. 

If one considers a ID long Josephson 0-7r-j unction 
(LJJ) made of alternating parts with positive and neg- 
ative critical currents (0 and 7r-parts), half-integer flux 
quanta (scmifluxonsi£) may spontaneously form at the 
boundaries between and 7r regions. 

Semifluxons are very interesting objects which are not 
yet studied in detail, first of all because up to now it 
was rather difficult to fabricate 0-7r-junctions. Recently 
several groups succeeded to demonstrate 0-7r-junctions 
based on various technologies: YBa2Cu307-Nb ramp 
zigzags^, grain boundary junctions based on tri- and 
tetra-crystals^ and Nb junctions based on an artificially 
created discontinuity^. Both and tt junctions from 
the Superconductor- Ferromagnet-Superconductor family 
were demonstrated by several groups&2*i2i but a 0-7T-LJJ 
was not reported yet. Semifluxons were observed using 
SQUID microscopy in different types of 0-7T-LJ Js2*ii*i2*i£. 

A single semifluxon, formed in a LJJ of length 
L 3> Aj with one 0-7r-boundary, is pinned at this 0-n- 
boundaryii*^ and can have positive or negative polarity 
carrying the flux +<J>o/2 or — $o/2, respectively. The bias 



current from up to ~I C (J c = j c wL is the "intrinsic" 
critical current, w is the junction width) cannot move 
the semifluxon but just changes its shaped. This prop- 
erty suggests to use semifluxons in information storage 
devices, classical or quantum. 

In the classical regime the polarity of a semifluxon 
(positive or negative) will encode a logical or 1. The in- 
formation encoding using semifluxons is somewhat more 
robust than using fluxons because semifluxons cannot 
"escape" as they are pinned at the discontinuity point. 
The switching between these states can be done by in- 
jecting a single flux quantum of proper polarity into the 
junction. 

In the quantum limit a semifluxon having two possible 
polarities is similar to a spin with two possible orienta- 
tions (up or down) . Therefore we often use spin notation 
to denote the polarity of the semifluxon, i.e., f or J,. 
Analogously, we denote the polarity of a fluxon as ff or 
JJ-. It seems that a semifluxon is an interesting candidate 
to realize a qubit. It still remains a challenging task to 
find out whether the semifluxon may stay in the super- 
position of both states and perform quantum tunneling 
between them or not. In comparison with a fluxon based 
qubit the one based on a semifluxon should be more ro- 
bust as both I and J, states represent the ground state 
of the system, while the fluxon of any polarity is an ex- 
cited state, the ground state being the one with constant 
phase. 

For both quantum and classical bits one needs a way 
to determine the final state of the semifluxon, i.e., read 
out its polarity. Imagine the simplest situation: a single 
semifluxon of unknown polarity in an annular LJ J2£. Let 
us inject a fluxon of a certain polarity into this Josephson 
ring somewhere far from the semifluxon position, e.g., us- 
ing a pair of current injectors^. If the polarities of the 
fluxon and of the semifluxon are different, the "annihi- 
lation" between fluxon and semifluxon will result in a 
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pinned semifluxon of the opposite polarity (we assume 
that the bias current is zero during fluxon injection). 
The I-V characteristic (IVC) of the resulting state has 
a rather large maximum supercurrent (depinning current 
of a semifluxon) 7J = 2/tt in normalized units (7 = I j 'I c ). 
On the other hand, if the fluxon and the semifluxon are 
of the same polarity, no "annihilation" takes place. If the 
bias current is applied, the fluxon starts moving passing 
through the semifluxon resulting in a finite voltage across 
the LJJ and rather low j c . In this case the depinning 
current 7JT is determined by the repulsion force between 
fluxon and semifluxon. The fluxon should overcome this 
pinning by the semifluxon to start moving around the 
ring. 

In this paper we propose a technique to test the po- 
larity of semifluxons by introducing test fluxons into the 
Joscphson ring. Measuring the current voltage character- 
istic we can determine the polarity of semifluxon(s) be- 
fore the fluxon was introduced. We also study resulting 
states which have different critical (depinning) currents 
and different dynamics. 

In section [n] we introduce the model which is used for 
numerical simulations. The numerical results arc pre- 
sented and discussed in section Hill Section llVl concludes 
this work. 



II. THE MODEL 

The dynamics of the Josephson phase in a LJJ consist- 
ing of alternating and 7r parts can be described by the 
ID perturbed sine-Gordon equation^ 

4>xx ~ 4>tt - sin0 = a(f> t - 7 - 6 xx (x), (1) 

where <f)(x, t) is the Josephson phase and subscripts x 
and t denote the derivatives with respect to coordinate 
x and time t. In Eq. Q the spatial coordinate is nor- 
malized to the Josephson penetration depth Aj and the 
time is normalized to the inverse plasma frequency w ~ ; 
a = is the dimcnsionless damping (/3 C is the 

McCumber-Stewart parameter); 7 = j/j c is the exter- 
nal bias current density normalized to the critical current 
density of the junction. The function 9(x) is a step func- 
tion which is 7r-discontinuous at all points where and tt 
parts join and is a constant equal to nn within each part 
(n is an integer). For example, 9(x) can be equal to zero 
along all 0-parts and 7r along all 7r-parts. 

It is clear from Eq. Q that 4>{x) is also 7r-discontinuous 
at the same points as 0(x). Therefore, we often call 
the points where and 7r parts join phase discontinuity 
points. 

Note that to describe 0-7r-LJJs other author a 2 ^^ 1 ? of- 
ten use directly the equation with alternating critical cur- 
rent density written for the continuous phase /j,(x, t) 

Hxx - P-tt ± shi/i = apit - 7(x). (2) 

The Eqs. Q and J3J) are, actually, equivalent and one 



can be obtained from the other by substitution (f>(x, t) = 
H(x,t) +8{x)L. 

In case of an annular LJJ one should use periodic 
boundary conditions (b.c.) to solve Eq. or In 
the case of a conventional annular LJJ without phase 
discontinuities the boundary conditions are expressed as 
follows 

<t>{L,t)=<t>{Q,t)±2im F , (3) 

where np is the number of flux quanta (Josephson vor- 
tices) trapped in the ring. Note that when there are no 
discontinuities \i = <fi, so that b.c. © holds for [i also. 

For the case of an annular LJJ with discontinuities 
the boundary conditions for <j) arc still given by Eq. J3J). 
This can be understood using the following gedanken ex- 
periment. Imagine that we start from the state with- 
out discontinuities <p{x) = 0. Then we slowly increase 
a discontinuity, e.g., by using a pair of closely located 
(5-injectorsi, at some point x = xo from the value to 
some value k. It is clear that the Josephson phase <j>{x) 
changes somehow on the length scale Xj in the vicinity 
of xq to compensate (to react on) this discontinuity, e.g., 
by forming a fractional vortex with the center at x = xq . 
In any case the phase (j>(x) and its derivative <p x (x) are 
smooth and continuous functions all along the junction 
except for the point x = Xq. Assuming that the disconti- 
nuity point does not coincide with x = or with x = L, 
we can write 4>(0) = 4>(L) and <f> x (Q) = 4> X (L). In the 
presence of additional fluxons trapped in the junction we 
get b.c. © even in the presence of the discontinuity 
points. 

The b.c. for fj, can be written, recalling that fi = <fi — 
6{x) 

H(L, t) = //(0,t) + ttosf + 2ttti f , (4) 

where «sf is the sum of all discontinuities (semifluxons) 
in the ring. Note that if 

N 

9(x) = KjHjx - Xj), (5) 

i=l 

where m is the i-th discontinuity and Ti(x) is a Heavyside 
step function, the expression for tisf is 

N 

nsF^^Ki, (6) 

i=i 

In this paper we will investigate the following two cases 

1. A fluxon is injected into the annular LJJ containing 
a negative semifluxon, i.e., np = +1, usf = — 1. 

2. A fluxon is injected into the annular LJJ containing 
a positive semifluxon, i.e., np = +1, Usf = +1. 
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III. NUMERICAL RESULTS 

The simulations were performed using StkJJ 
software^ and were confirmed by an independently writ- 
ten program 2 ^. 

We use a = 0.1 for all results reported here. This 
value is not very high, so that it allows to observe some 
dynamics. On the other hand, it is not very low as the 
majority of 0-7r-j unctions has rather high damping. In 
the case a > 1 the static results related to the reading 
out the state of a semifluxon hold, but no dynamical 
effects such as fluxon steps can be observed. 

In order to visualize and to understand the fluxon and 
semifluxon dynamics we plot their trajectories on the 
(x, t) plane. Usually, to track the trajectory of a fluxon 
one tracks the trajectory of its center. Since at a given in- 
stant of time the phase /z(x), corresponding to the fluxon 
solution, changes from at x — > — oo to 2n at x — > +oo, 
it is assumed that the center of a fluxon coincides with 
the point where the phase fi = n. Since the phase is de- 
fined modulo 2tt, in the general case the center of fluxons 
is situated at points where (i = n + 2nk with integer k. 
To distinguish between fluxons and antifluxons, we also 
check the sign of the phase derivative fx x (magnetic field) 
at the point where fj, = tt + 2i:k. If the sign is positive, 
then it is a fluxon and we plot its position as a black 
point on x-t plane. If \i x < 0, then we plot it as a gray 
points on the x-t plane. 

When we deal with semifluxons, the idea is the same, 
but since the phase of a semifluxon changes from to 7r 
( mod 7r) we have to define the centers of semifluxons 
as the points where fi = ^ + irk. The sign of fi x at the 
points where /i = ^ + irk is used to distinguish between 
semifluxons of positive and negative polarity. With this 
definition, every 27r-fluxon carrying the integer flux $o 
results in two points on the x-t plane, i.e., its trajectory 
will be represented by a double line. 

Below we present numerical results obtained for an an- 
nular junction of length L = 8Xj. 

A. One semifluxon in a ring 

First we investigate numerically an annular LJJ with 
only one phase discontinuity point. 

Here and below, without loosing generality, wc assume 
that the injected fluxons have positive polarity ff [np = 
+1 in gj}]. 

1. Negative semifluxon 

If initially the semifluxon has negative polarity j 
["■sf = —1 hi Q] the fluxon is attracted by the semi- 
fluxon and, in the absence of a bias current, they "anni- 
hilate", resulting in a positive semifluxon f. 

In Fig. Ola) one can sec the "annihilation" process of 
the semifluxon situated at x = 3 and the fluxon injected 



at x = 6 (it corresponds to two lines at x 5 and x « 7 
at t = 0) for zero bias current 7 = 0. If we trace the IVC 
of the state after annihilation, wc get the curve shown 
in Fig. |21 by solid black symbols. As a horizontal axis 
we use the fluxon velocity u normalized to the Swihart 
velocity and proportional to the voltage across the junc- 
tion. In this notation a single fluxon has an asymptotic 
velocity 1, two fluxons have u — > 2 as 7 grows, etc.. At 
u = the state of the system is the semifluxon f with 
maximum supercurrent (semifluxon depinning current) 
jl = 2/tt w 0.63 14 i 20 i 21 i 22 . When the bias current ex- 
ceeds this value the system switches to the McCumber 
branch. By sweeping the bias current back one may trace 
the step with asymptotic velocity u = 3, corresponding 
to the state 'fMNI-J,, which is formed probably due to the 
topological instability of the solution similar to the for- 
mation of zero- field steps in conventional LJJ. The tra- 
jectories corresponding to the bias point d can be seen in 

Fig.ntd). 

What happens if we inject a fluxon while having a non- 
zero bias current? If the bias current as well as the dis- 
tance between fluxon injection point and semifluxon are 
large enough, the fluxon approaches the semifluxon with 
rather high velocity and annihilation will not take place 
- the fluxon will simply pass through the semifluxon. 
Thus, injecting a fluxon at a finite bias current we can 
trace an IVC which is shown in Fig. [3 by open symbols. 
The semifluxon trajectories at the bias points b and c are 
shown in Fig.^b)-(c), accordingly. One can see that in 
average a semifluxon is shifted by the bias current away 
from the discontinuity point and oscillates around this 
new equilibrium position as the fluxon bumps it. Com- 
paring Fig. ^b) and (c), one can also notice that the 
fluxon's double line is more tight in (c) which is a result 
of relativistic contraction. Note that after the bias cur- 
rent is reduced to zero the annihilation takes place any- 
way and the system returns to the IVC corresponding 
to the state |, without a moving fluxon, so that further 
sweeping of 7 shows only the curve drawn by solid black 
symbols in Fig. |21 

Actually one can also trace the second fluxon step with 
asymptotic velocity u = 2. It corresponds to the state 
f|"JJ.f and is shown in Fig. |21 by solid gray symbols. We 
were able to find this state only starting from the point 
e at 7 = 0.3 with a fluxon situated at x = 5 and an 
antifluxon at x = 8. It is impossible to visualize this 
mode just by sweeping the bias current since this step is 
shadowed by t, fUT-IJ-l and -ft J, steps. 



2. Positive semifluxon 

Now we consider a semifluxon of a positive polarity ini- 
tially present in the LJJ. The injected fluxon and semi- 
fluxon repel each other and no annihilation can occur. 
Here we use the boundary conditions gj with hp = 1 
and /csf = 1 to obtain the IVC. By applying a small bias 
current we push the fluxon along the ring so it approaches 



4 





mm 

(a): 

) ; 




A i 1 " 



2 4 6 8 
coordinate x 




coordinate x 



coordinate x 



1 1 I 1 i 1 

2 4 6 
coordinate x 





■< i 1 i 1 

(e). 

i 4 


y 




A 

> i 1 


1 1 1 1 1 



8 2 4 6 
coordinate x 



FIG. 1: Trajectories of fluxons (double line) and semifluxons (single line) corresponding to annihilation ^ + | = |at7 = (a); 
dynamics in the state ftj at the first fluxon step at 7 = 0.3 (b); at 7 = 0.6 (c); dynamics at the third fluxon step (state fnHH) 
at 7 = 0.6 (d); dynamics at the second fluxon step (state 1HI-T) at 7 = 0.3. The vertical dashed line shows the position of the 
discontinuity. The vertical dotted line shows the position of the center of the static semifluxon at the corresponding bias. In 
(a) for 7 = both lines coincide. 
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FIG. 2: IVC after injection of a fluxon into a 0-7T-LJJ con- 
taining a semifluxon of negative polarity. The solid black 
symbols show the IVC after annihilation, open symbols show 
the IVC which can be traced if injection takes place at finite 
2/tt > 7 > y r ~ 0.21. Dotted line shows the position of the 
McCumber branch (uniform phase- whirling state). 



FIG. 3: Current voltage characteristic after injection of 
fluxon into 0-7T-LJJ containing a semifluxon of positive po- 
larity. All (even and odd) fluxon steps can be traced. The 
maximum supercurrent 7 max ~ 0.21. The trajectories at the 
bias points a-d are shown in Fig. Bla)-(d). 



the semifluxon. For a small value of the bias current the 
situation is static as the driving force of the bias cur- 
rent can be compensated by the repulsion force between 
the fluxon and the semifluxon. By increasing the bias 
current, the fluxon moves closer to the semifluxon and 
at some critical value of the bias current overcomes the 
maximum possible repulsion force and passes through the 
semifluxon. After this, the fluxon keeps moving around 
the ring, bumping the semifluxon once per cycle. The 
current voltage characteristic of this state is shown in 
Fig. 13 One can see that the maximum supercurrent 
7^ = ss 0.21 (see Appendix |A"|) . corresponding to the 
maximum possible repulsion force between the fluxon and 
the semifluxon, is considerably smaller than 7J = 2/tt of 
a single semifluxon. Thus, the two situations can be dis- 
tinguished. 



The trajectories corresponding to the dynamics of the 
system are shown in Fig. 0] for several bias points marked 
in Fig. |21 We see that at the first fluxon step, see 
Fig. 01b), the fluxon (double line) moves progressively 
colliding with the semifluxon. The semifluxon corre- 
sponds to a more or less vertical line shifted to the right 
from the discontinuity point by the bias current. Note, 
that in Fig. 0fa) we can see essentially the same dynam- 
ics, but the semifluxon is much more delocalized, prob- 
ably because bias point a corresponds to the resonance 
which can be seen on the IVC in Fig. [S] We believe 
that this resonance may be related to the eigen-modes 
of the semifluxon and to the Cherenkov emission of the 
plasma waves which are excited when the fluxon period- 
ically bumps the semifluxon. The emitted plasma wave 
forms a standing wave which interacts with a fluxon and 
results in the resonance on the IVC. 
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FIG. 4: Trajectories of semifluxons in the state -f|-"f (first fluxon step) for 7 = 0.2(a) and 7 = 0.6 (b) and in the state •ft-frj. for 
7 = 0.3 (c) and 7 = 0.6 (d). The vertical dashed line shows the position of the discontinuity. The vertical dotted line shows 
the position of the center of the static semifluxon at the corresponding bias. In (a) for 7 = both lines coincide. 



Thus, the polarity of the semifluxon can be probed by 
inserting a test fluxon of known polarity and measuring 
the IVC. If a fluxon and a semifluxon have the same 
polarity one should expect I max = ^ « 0.21/ c and the 
appearance of the first fluxon step. If a fluxon and a 
semifluxon have different polarities, one should expect 
I max = — I c m 0.63/ c and no first fluxon step provided 
the injection was made at 7 = 0. 

Another difference between these two configurations is 
the value of the retrapping current y r . As one can see 
from Figs. |2 and |3] the value of 7 r in the ffl = T state is 
at least twice larger than for the ^ state. 

An interesting observation can be made, if we trace all 
fluxon branches of the IVC. In a usual annular LJJ with 
only one trapped fluxon (without a semifluxon) one can 
observe fluxon steps at V = nV\ with only odd numbers 
n: the first one appears due to the motion of a single 
fluxon, the third one corresponds to an additional fluxon- 
antifluxon pair generated, the fifth one to two fluxon- 
antifluxon pairs, etc.. 

In an annular LJJ containing a fluxon and a semifluxon 
of the same polarity the situation is different. One can 
observe fluxon steps with any integer n. This is especially 
easy to observe for the ■frf state. The reason for this 
becomes clear after analysis of the trajectories shown in 
Fig. 21 for several bias points marked in Fig. [21 

We see that at the second step, there are two positive 
fluxons moving in the same direction and the semifluxon 
became of negative polarity ififj,. This means that the 
semifluxon flipped, and has emitted one positive fluxon 
— the process opposite to annihilation. Now, since there 
are two fluxons which can move and one negative semi- 
fluxon which cannot, the asymptotic voltage of the step 
is equal to 2V\. Still, the system can generate fluxon- 
antifluxon pairs which, together with semifluxon flipping, 
will result in fluxon steps for all integer n. 

Actually, a similar effect can be observed when initially 
the semifluxon had negative polarity. Then starting from 
the state | one could trace the fluxon steps corresponding 
to the states 1U, 1T1TJJ4 (sec Fig.|2J). The trajectories 



for some of the bias points can be seen in Fig. 2] 

We also note that when the semifluxon polarity 
changes, the shift of the average position of the semi- 
fluxon changes too. This can be visualized using low 
temperature scanning electron microscopj*2£. 

Analyzing the trajectories while the bias point moves 
along the second fluxon step we discovered that the two 
fluxons are moving more or less equidistantly at a low 
value of the bias current, i.e., at the bottom of the step, 
as can be seen in Fig. 0{c). When the bias current 
increases up to 7 « 0.45 . . . 0.50 or above, the fluxons 
bunch together as shown in Fig. 0Jd) . The bunching ap- 
pears because the plasma waves, emitted during fluxon- 
scmifluxon collisions, result in an effective attracting po- 
tential between two fluxons. We do not discuss this phe- 
nomenon in detail here, but note that similar bunching 
was observed in JJ arrays^ or in stacks of LJJgS. In our 
case, the role of periodic obstacle is played by a semi- 
fluxon. Generally speaking, the plasma wave emission 
has a Chcrcnkov origin and is simply related to the pe- 
culiar dispersion relation for plasma waves in the system 
under question^. 



B. Two semifluxons in a ring 

In the case of more than one semifluxon of unknown 
polarity, the above read out procedure cannot give infor- 
mation on the polarity of each of the semifluxons, but 
allows to distinguish between states (a) ||, (b) J. J., and 
(c) |J, or J.1". In other words, we can find out the total 
flux hold by semifluxons. Injecting the first fluxon in a 
similar fashion as above, we measure the IVC. If the crit- 
ical current is low, the fluxon does not annihilate with 
semifluxons and keeps moving around the ring. This can 
happen only if in the initial state both semifluxons had 
the same polarity as the injected fluxon. If, after in- 
jection of a fluxon, I c is still high, this means that the 
injected fluxon has flipped one of the semifluxons. In this 
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case we inject a second fluxon and so on. If I c becomes 
lower after the first injection, the initial state has been 
It; if I c becomes lower after second injection, the initial 
state has been fj or J, j; an d if I c becomes lower after the 
third injection, the initial state has been ||. 



IV. CONCLUSIONS 

We have shown that by introducing a test fluxon of 
known polarity into the LJJ with a semifluxon of un- 
known polarity we can destructively read out the semi- 
fluxon state. In principle, the read out can be made 
non-destructive, if after reading out the semifluxon state 
we introduce another fluxon of opposite polarity so that 
the system returns to its initial state (as before read out) 
as a result of fluxon-fluxon or fluxon-semifluxon annihila- 
tion. In case of two or more semifluxons one can read out 
the total number of positive and negative semifluxons. 
This technique is also applicable to arbitrary fractional 
K-vortices, but not very close to k = or k = 2ix where 
depinning currents for both N — and N = 1 approach 
zero. 

We have also investigated fluxon steps in an annular 
LJJ with a 7r-discontinuity. We found that one can trace 
the fluxon steps corresponding to all integer n, instead 
of only odd n like in conventional LJJ with one trapped 
fluxon. We have observed a smooth transition to the 
state of two bunched fluxons in the f|-f|-J, state. 

Further, in appendix lAl we have derived the Eq. (|A9|) 
which gives the depinning current of N fluxons pinned 
by an arbitrary k- vortex. We discovered that the biggest 
obstacle for a fluxon is a fractional Ki-vortex with n\ w 
0.86I7T rather than a semifluxon. The formula (|A10|) 
allows to compute the size kjv of the fractional vortex 
which is the biggest obstacle for N fluxons trying to pass 
it. 
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APPENDIX A: DEPINNING OF FLUXONS BY 
AN ARBITRARY FRACTIONAL VORTEX 



where 



6{x) 



for x < 0, 
— k for x > 



We write separately the Eqs. for the part of the junc- 
tion to the left and to the right from discontinuity, sit- 
uated at x = 0. After integration one arrives to the 
following two equations. 



2/i(/0 



±2^/C^) 



cos n — 7/i 
±2VC 2 (7) -cos(/x 



x < 0; (A2a) 
k) — 7//, x(A2b) 



Assuming that 

Ma(±oo) 
n(-oo) 
/i(+oo) 



0; (A3) 
arcsin 7; (A4) 
arcsin 7 + 2-kN + k, (A5) 



we arrive to the following expressions for C± and C2' 



Ci(7) = yl — 7 2 + 7 arcsin 7; 



(A6) 



C 2 {l) = y/l - 7 2 + 7 arcsin 7 + (2irN + k)j, (A7) 

We use the phase plane analysis^ to find possible 
static configurations. The trajectories on the phase plane 
Ui,2{^) corresponding to Eqs. ljA2() are shown in Fig. 
for 7 = 0, 7 < 7 C and 7 = j c . At 7 = the fluxon 
and fractional vortex are separated by a large distance. 
The fluxon corresponds to the trajectory between point 
L (x = —00) and point M. The fractional vortex cor- 
responds to the trajectory between point M and point 
R (x = +00) and contains a point P where the black 
and gray trajectories intersect, i.e., the phase crosses the 
0-7r-boundary. From Fig. Efa)-(c) one can see that with 
increasing 7 the intersection point P shifts until, at the 
critical value of 7 = j c the trajectories just touch each 
other at point P, as shown in Fig. EJc). For larger 7, no 
intersection is possible and the static solution does not 
exist. 

The critical value of 7 at which the switching between 
trajectories is still possible is defined by the following 
conditions: 



2/1 Cu) = y2(fj.);y' 1 {fj.) = y^n), 



(A8) 



This conditions are satisfied for [i = ^(37r + k), which 
leads us to the final result: 



Since one can study experimentally arbitrary n- 
vorticesi, here we calculate the depinning current for a 
chain of N fluxons that are pinned by and are trying to 
pass an arbitrary K-vortex in a LJJ of infinite length. 

The static version of Eq. J2J for arbitrary discontinuity 
k is: 



Hxx = sin[^ - 6(x)] - 7, 



(Al) 



1c(k) 



sin — . 

2ttN + k 2 



(A9) 



The plot of this dependence for different N is shown 
in Fig. The j c (k) has reasonable limiting behavior 
(absence of pinning for N > 0) for k — > and n — > 
2t:. For depinning of a fluxon (N — 1) by a semifluxon, 
7c(tt) = 2/3-7T w 0.21. It is interesting that 7 C («) has a 
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Note, that in the limit N — > oo, kjv - * T. 

Our result 1A9|I is also valid in the case N = and gives 
the "depinning" current of the K- vortex itself. The result 
(|A9|) coincides with the expression obtained recently for 
the critical current of an annular junction in the presence 
of a /{-discontinuity of the phase (current dipole) created 
by injectors^. In fact one can even say that our result 
1.0 



c 

0) 

o 

O) 

c 
'c 
c 

Q. 

CD 

Q 




0.5 1.0 1.5 

fractional vortex k/ti 



2.0 



FIG. 6: The dependence 7 C («) calculated for different N using 
Eq. <A9t . The dashed line shows k, = ir corresponding to the 
semifluxon. Three dotted lines show the values Ki, ft2 and K3 
which correspond to the fractional vortices which cause the 
maximum possible pinning for N = 1, 2, 3 fluxons. 



FIG. 5: The trajectories on the phase plane corresponding 
to Eqs. ljA2fl . Black color corresponds to yi (O-part), the 
gray color corresponds to j/2 (7r-part). The trajectories are 
presented for a fluxon pinned by a semifluxon (N = 1, k = 7r) 
and (a) 7 = 0, (b) 7 = —0.1 and (c) 7 = 0.21. Arrows 
are shown to guide the eye. Arrows show the path which 
corresponds to going from x — —00 (point L) to x = +00 
(point R). 



maximum not at k = n but a bit shifted, at K\ = 0.8617T, 
i.e., a Ki-vortex is the biggest obstacle for a fluxon. For 
two and three fluxons the biggest obstacles are vortices 
with K2 = 0.918-7T and K3 = 0.9427T. All kn can be found 
from the equation 



tan — = h Ntt. 

2 2 



(A10) 



(|A9(1 and the one of Ref. [2(j are the same at least if one 
considers annular LJJ. In fact, if we denote k' = 2irN + k, 
we can rewrite (IA9I) as 



7cOO = 



2 sin ^ 



(All) 



exactly as in Ref. |2fJ. Physically this means that when 
one creates a large discontinuity k! = 2ttN+k in an annu- 
lar LJJ, it automatically relaxes into N fluxons trying to 
pass through the fractional K-vortex under the action of 
bias current. The interesting result is that the Eq. (|A9|) 
for an infinitely long linear junction and Eq. I|A11|) for 
an annular one coincide regardless of "fluxon crowding" 
which may take place in annular LJJ because of its finite 
length. 
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We note that such annular junctions cannot be constructed 
using natural or tt junctions because one part cannot be 
simultaneously and n. There should be always an even 
number of 0-7r-joints. In an experiment the proposed sit- 
uation can be realized using artificial 0-7T-LJJ based on 
conventional superconductors with a pair of 5-injectors-. 
This case of only one discontinuity in a ring is very inter- 
esting as it allows to go beyond the possibilities offered by 
nature. 
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